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Abstract 

A modified version of the M0ller-Plesset approach for obtaining the correlation 
energy associated to a Hartree-Fock ground state is proposed. The method is tested 
in a model of interacting fermions that allows for an exact solution. Using up to 
third order terms improved results are obtained, even in the limit of loosely bound 
particles. 



The study of molecules and larger systems is seriously constrained by their 
many body nature. Several approximation schemes have been devised over 
the years, among which the Hartree-Fock (HF) method is one of the oldest 
and most fruitful. Because it treats interactions in a mean field way particle 
correlations are left out, however, a shortcoming that can limit severely the 
validity of its results. One may improve over HF by treating correlations as a 
perturbation. In the so-called M0ller-Plesset method, a Rayleigh-Schrodinger 
perturbative expansion that is naturally suggested by the same structure of 
the HF solution[l,2,3,4] is adopted. This formalism reduces the correlation 
energy to an infinite series in the perturbation, of which only the first few 
terms need to be computed in practice. This scheme has been used for a long 
time as a good starting point to study correlation effects in molecular systems. 
However, this method is useful only if the perturbation series is rapidly conver- 
gent, which is not always true [5]. Failures of the M0lIer-Plesset method have 
been documented even for small molecules [5,6,7]. In this paper we present a 
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variation of the M0ller-Plesset approach that appears to give accurate results 
in low order perturbation schemes even when particles are loosely bound, such 
as in chemical bonds. 



Consider the many-particles Hamiltonian 

w = EM0 + E v (m)» (i) 

i i,j>i 

where h(i) = hfa) is the sum of the one-particle kinetic energy plus external 
potential energy, and v(i,j) = vfa, Tj) is the two-particles interaction energy. 
The HF approximation to the ground state of a system of identical fermions 
leads to a variational wave function <E> in the form of a Slater determinant, 
constructed with one-particle orbitals that satisfy the equations 

h(i) + E " M*))] M*) = £n0n(i), (2) 

b 

where Tb and /C& are the Coulomb and exchange operators, respectively. Here 
and in what follows summation over indices a, b, c run over all occupied states, 
while n runs over all possible states. The eigenvalues satisfy the relation 

e n = (n\h\n) + J2(nb\\nb), (3) 

b 

where 

(mn\\mn) = (mn\mn) — (mn\nm), (4) 
(mn\pq) = J m (l)*0 n (2)^(l, 2)^(1)^(2)^2. (5) 

In the above expression 1 and 2 represent the one-electron variables of coor- 
dinate and spin. The energy of the HF state is 



a Z a,b 

= £e a - V(a%6>, (6) 

a Z a ,b 

where to obtain (6) we have used Eq. (3). The operator 

7-Chf = E e nclc n - i ^2(ab\\ab) (7) 

n ^ a,b 

is then a natural choice for an effective HF hamiltonian, diagonal in the HF 
orbitals and having as ground state the HF energy E HF . In Eq. (7) (c n ) is the 
creation (annihilation) operator of a particle in the state <f> n . This hamiltonian 
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operator is the starting point to construct a perturbation expansion for the 
correlation energy using as perturbation Vhf — 'H — 7~Lhf [4]. 



A second choice for a hamiltonian is also possible, however. Equation (6) may 
be written in the form, 

^F = £(e„-^ (8) 



2 b 



suggesting as an alternate hamiltonian 



Umhf = £(e„ - \ Y,( nb \\ nb ))c{c n , (9) 

n A b 

also diagonal in the HF orbitals. The correlation energy can again be described 
in terms of a perturbation, this time of the form Vmhf = 7~L — 7~Lmhf- Note 
that while Eq. (7) involves an overall constant, Eq. (9) substracts from each 
single particle energy a different correction. Although both forms yield the 
same ground state energy, the energy of excited states are different, affecting 
the various orders contributions in a perturbation expansion. For instance, the 
second order correction to the HF ground state energy has the form [4] 

$/ &HF — -CV 4 abrg £j HF — h,& 

where in the numerator we have used the fact that the excited state $' is or- 
thogonal to the ground state. This is the first finite correction in the Rayleigh- 
Schrodinger perturbation expansion since for either choice of hamiltonian the 
first order term vanishes. Only doubly excited states $' = $^ (<j) a replaced by 
<ft r and <ftb replaced by <p s in the Slater determinant $, with e r , e s above the 
Fermi energy) contribute [4]. The choice of hamiltonian affects the excitation 
energies in the denominator, which have the form 

Ehf — E$rs = ^ , (11) 

£ sj 7~Lo = Hmhf 

where e n = e n — \J2b{ n H nb ) = |( e « + (n\h\n)). Replacing in Eq. (10) both 
forms clearly lead to different results. A similar analysis of the third and higher 
order terms in the perturbation expansion yields again different results. For 
example, in the third order correction, besides the change of denominators in 
the standard M0ller-Plesset expression [4], an additional term AE^ appears, 
which can be cast in the compact form, useful in computations, 

A£( 3 ) = -E^ 

1 x ^ {hrr + h ss - h aa - h bb )\(ab\\rs)\ 2 

*L {ir + is-ia-h) 2 ' 1 j 
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where h nn = (n\h\n). The first element of Eq. (12) cancels the second order 
energy correction. However, the second element of Eq. (12) is also second order 
in the interaction. Replacing Eqs. (10) and (11) in (12) and performing simple 
operations AE^ can be put in a form that reveals explicitly the third order 
character of the correction. 

Table 1 

Ground state energies of a system of two harmonically confined spin 1/2 particles 
interacting through a harmonic potential of strength k. Several approximations are 
included: Hartree-Fock (HF), Moller-Plesset perturbation theory of orders n = 2,3 
(MPn), modified MPn (MMPra), as well as the exact values. Repulsive interactions 
are represented by negative values of the elastic constant k. 
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k 


1 IT" 

HI 


TV r i ~> ■ i 

MP2 


MP3 


TV it TV TTiO 

MMP2 


TV TTV r 1 ~> O 

MMP3 


Exact 


-0.25 


1.732 


1.655 


1.836 


1.702 


1.710 


1.707 


-0.24 


1.744 


1.681 


1.803 


1.717 


1.724 


1.721 


-0.22 


1.766 


1.725 


1.784 


1.745 


1.750 


1.748 


-0.20 


1.789 


1.760 


1.791 


1.772 


1.776 


1.775 


-0.18 


1.811 


1.791 


1.808 


1.798 


1.801 


1.800 


-0.16 


1.833 


1.819 


1.828 


1.823 


1.825 


1.825 


-0.09 


1.908 


1.905 


1.906 


1.905 


1.906 


1.906 


-0.04 


1.960 


1.959 


1.959 


1.959 


1.959 


1.959 


-0.01 


1.990 


1.990 


1.990 


1.990 


1.990 


1.990 


0.00 


2.000 


2.000 


2.000 


2.000 


2.000 


2.000 


0.04 


2.040 


2.039 


2.039 


2.039 


2.039 


2.039 


0.16 


2.154 


2.150 


2.149 


2.149 


2.149 


2.149 


0.36 


2.332 


2.319 


2.314 


2.316 


2.313 


2.311 


0.64 


2.561 


2.534 


2.522 


2.525 


2.516 


2.510 


1.00 


2.829 


2.784 


2.762 


2.767 


2.749 


2.732 



In order to assess the convenience of either formulation for obtaining correc- 
tions due to correlations we have solved a system of two spin-1/2 particles in a 
harmonic potential, interacting through a harmonic force. The reason for the 
choice is that this is an interacting system involving identical fermions that 
may be solved exactly [8]. The Hamiltonian is 

H = \ (-A x + r?) + \ (-A 2 + rl) + - r 2 ) 2 , (13) 
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where the coordinates and the energy are given in the oscillator units of the 
confinement potential. The exact ground state energy is 

£ =1 + VT+2k. (14) 

The model contains the parameter k that allows the study of attractive (k > 0) 
as well as repulsive (k < 0) interactions. Equation (14) shows that k = —0.5 
is the lowest value for which a bound state exists. 

We have solved the problem in two dimensions, using a basis of noninteracting 
harmonic oscillators eigenfunctions 

( l>n x ,n y {Xiy)=VnS X ) i Pny{y), < U x + U y < 5, (15) 

where the tp(x) are the usual one-dimensional harmonic oscillator eigenfunc- 
tions. The set (15) is an exact solution when k = 0. The HF solution is obtained 
by solving the self-consistent-field equation for closed shell configurations^]. 
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Fig. 1. Correlation energy for a system of two harmonically confined spin 1/2 par- 
ticles interacting through a harmonic potential of strength k. The exact result is 
included, as well as corrections up to second and third order for the M0ller-Plesset 
(MP) and modified M0ller-Plesset (MMP) choices of zeroth-order hamiltonian. 

Table 1 shows the ground state energy for the Hamiltonian (13) calculated us- 
ing the standard M0ller-Plesset perturbation theory of order n = 2, 3 (MPn), 
and our modified form, MMPn. Figure 1 shows the exact correlation energy 
E corr = Ehf — E for different values of k, together with results obtained for 
the two choices of HF hamiltonian, in second and third order of perturbation 
theory. Notice that MMPn yields better results throughout. Notice also that 
the usual MPn fails badly in both orders of approximation when the system 
becomes more loosely bound, as k approaches the critical value -0.5. In fact, 
convergence problems prevent solving the HF self-consistent equations for k 
beyond -0.25. By contrast, MMPn continue to be good approximations even 
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in this range of k. This strongly suggests that the modified perturbation the- 
ory may be more suitable in treating outer shells of bound systems, such as 
electrons participating in chemical bonds. Investigation of this ansatz is in 
progress. 

A well known feature of HF theory is that the many body ground state energy 
is not the bare sum of energies of filled single particle orbitals. Interactions 
are counted twice and this overestimation is corrected for by substracting 
the constant explicit in Eq. (3). The remarkable improvement obtained over 
the usual M0ller-Plesset approach in our test example may be traced to the 
fact that this latter method does not correct the single particle energies for 
such effect. In fact, the energy denominator appearing in the perturbative 
corrections to all orders in such case may grossly depart from the true two- 
particle excitation energy that the numerator in the expression is weighting. By 
contrast, in the method proposed here each self energy is corrected accordingly. 
It is hoped that our results will stimulate the use of the proposed method 
in situations where corrections to the Hartree-Fock approximations may be 
necessary. 
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